Planets in Stellar Clusters Extensive Search. III. A search for 
transiting planets in the metal-rich open cluster NGC 6791. 1 



in 

o 
o 

(N 
o 

Q 

> 

in 



o 
in 
o 

S3 

Pl,: 

6 

CO 

C3 



X 



B. J. Mochejska 2 

Purdue University, Department of Physics, 525 Northwestern Ave., West Lafayette, IN 47907 

bmoche j s@cf a . harvard . edu 

K. Z. Stanek, D. D. Sasselov, A. H. Szentgyorgyi, G. A. Bakos 2 ' 3 , J. Devor, V. Hradecky, 

D. P. Marrone, J. N. Winn 2 & M. Zaldarriaga 
Harvard- Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 
kstanek, sasselov, saint, gbakos, j devor, vhradecky, dmarrone , jwinn, 
mzaldarriagaOcf a. harvard. edu 



ABSTRACT 

We have undertaken a long-term project, Planets in Stellar Clusters Extensive Search 
(PISCES), to search for transiting planets in open clusters. In this paper we present the re- 
sults for NGC 6791 - a very old, populous, metal rich cluster. We have monitored the cluster 
for over 300 hours, spread over 84 nights. We have not detected any good transiting planet 
candidates. Given the photometric precision and temporal coverage of our observations, and 
current best estimates for the frequency and radii of short-period planets, the expected number 
of detectable transiting planets in our sample is 1.5. We have discovered 14 new variable stars in 
the cluster, most of which are eclipsing binaries, and present high precision light curves, spanning 
two years, for these new variables and also the previously known variables. 
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We have undertaken a long-term project, Plan- 
ets in Stellar Clusters Extensive Search (PISCES), 
to search for transiting planets in open clusters. 
To date we have published a feasibility study based 
on one season of data for NGC 6791 (Mochejska 
et al. 2002, hereafter Paper I). We have also pre- 
sented a variable star catalog in our second target, 
NGC 2158, based on the data from the first ob- 
serving season (Mochejska et al. 2004, hereafter 
Paper II). 

In this paper we present the results of a 
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search for transiting planets in the open cluster 
NGC 6791 [(a, 6)2000 = (19 fe 20.8 m , +37°51')]- It 
is a very populous (Kaluzny & Udalski 1992), very 
old (r=8 Gyr), extremely metal rich ([Fe/H]=+0.4) 
cluster, located at a distance modulus of (m-M)y 
= 13.42 (Chaboyer, Green & Liebert 1999). 

Stars hosting planets are known to be, on the 
average, significantly more metal rich than those 
without (Santos et al. 2001, 2004). Two scenarios 
have been proposed to explain this phenomenon. 
Some studies favor a "primordial" metallicity en- 
hancement, i.e. reflecting the original metallicity 
of the gas from which the star formed (Santos et 
al. 2004; Pinsonneault et al. 2001). In this sce- 
nario planet formation would be more prolific in 
a metal-rich environment (Ida & Lin 2004). Oth- 
ers suggest that the host stars were enriched by 
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the infall of other giant gas planets (Lin 1997) or 
small planetary bodies like asteroids (Murray & 
Chaboyer 2002). 

The observed lack of planets in the core 
(Gilliland et al. 2000) and the uncrowded outer re- 
gions (Weldrake et al. 2005) of the low mctallicity 
([Fe/H]= -0.7) globular cluster 47 Tuc suggests 
that the source of the metallicity enhancement in 
planet hosts is most likely "primordial". Open 
clusters offer the possibility of observing a large 
number of stars with the same, known a priori 
metallicity. NGC 6791, with its high metallicity 
and large number of stars, seems particularly well 
suited as a target for transiting planet search. 

Targeting open clusters also eliminates the 
problem of false detections due to blended eclips- 
ing binary stars, which are a significant contam- 
inant in the Galactic field searches (over 90% of 
all candidates; Konacki et al. 2003; Udalski et al. 
2002a, 2002b). Blending causes a large decrease 
of the depth of the eclipses and mimics the transit 
of a much smaller object, such as a planet. As op- 
posed to dense star fields in the disk of our Galaxy, 
open clusters located away from the galactic plane 
are sparse enough for blending to be negligible. 

There are two key elements in a survey for tran- 
siting planets. The most commonly emphasized 
requirement is the high photometric precision, at 
the 1% level. The more often overlooked factor is 
the need for very extensive temporal coverage. 

Extensive temporal coverage is important be- 
cause even for planets with periods between 1 and 
2 days, the fractional transit length is only ~5% 
of the period, and it drops to ~2% for periods 2- 
10 days. During the remaining 95-98% of the pe- 
riod the system is photometrically indistinguish- 
able from stars without transiting planets. To 
our best knowledge, PISCES is the most exten- 
sive search for transiting planets in open clusters 
in terms of temporal coverage with aim tele- 
scope. 

NGC 6791 has been previously searched for 
transiting planets by Bruntt et al. (2003), who 
found three transit-like events and seven other 
lower probability events which may possibly be 
due to instrumental effects. Of the three best can- 
didates, none exhibited more than one transit and 
only one is located on the cluster main sequence. 
Bruntt et al. (2003) used the 2.5 m NOT telescope, 



which allowed them to obtain higher photomet- 
ric precision and denser time sampling, but their 
temporal coverage was much inferior to ours: ^24 
hours spread over 7 nights, compared to our > 300 
hours, collected over 84 nights. 

The paper is arranged as follows: §2 describes 
the observations, §3 summarizes the reduction 
procedure, §4 outlines the search strategy for tran- 
siting planets, §5 gives an estimate of the ex- 
pected number of transiting planet detections, §6 
describes the candidates previously reported by 
Bruntt et al. (2003) and §7 contains the variable 
star catalog. Concluding remarks are found in §8. 

2. Observations 

The data analyzed in this paper were ob- 
tained at the Fred Lawrence Whipple Observa- 
tory (FLWO) 1.2 m telescope using the 4Shooter 
CCD mosaic with four thinned, back side illumi- 
nated AR coated Loral 2048 2 CCDs (Szentgyor- 
gyi et al. in preparation). The camera, with a 
pixel scale of 0'.'33 pixel -1 , gives a field of view 
of 11. '4 x 1 1 '4 for each chip. The cluster was 
centered on Chip 3 (Fig. 1). The data were 
collected during 84 nights, from 2001 July 9 to 
2003 July 10. A total of 1118 x 900 s R and 
233 x 450 s T^-band exposures were obtained. The 
F-band dataset was supplemented with 93 x 450 
s exposures collected between 19 September 1998 
and 5 November 1999 (previously analyzed by 
Mochejska, Stanek & Kaluzny 2003). 

3. Data Reduction 

3.1. Image Subtraction Photometry 

The preliminary processing of the CCD frames 
was performed with the standard routines in the 
IRAF cedproc package. 1 

Photometry was extracted using the ISIS image 
subtraction package (Alard & Lupton 1998; Alard 
2000), as described in detail in Paper I. 

The ISIS reduction procedure consists of the 
following steps: (1) transformation of all frames 
to a common (x, y) coordinate grid; (2) construc- 
tion of a reference image from several of the best 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which arc operated by the Association of 
Universities for Research in Astronomy, Inc., under coop- 
erative agreement with the NSF. 
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Fig. 1 — Digital Sky Survey image of NGC 6791 
showing the field of view of the 4Shooter. The 
chips are numbered clockwise from 1 to 4 starting 
from the bottom left chip. NGC 6791 is centered 
on Chip 3. North is up and east is to the left. 



Fig. 2. — V/V — R color-magnitude diagram for 
Chip 3, centered on NGC 6791. Newly discov- 
ered eclipsing binaries are plotted with circles, 
other periodic variables with squares and the non- 
periodic variable with a triangle. 



exposures; (3) subtraction of each frame from the 
reference image; (4) selection of stars to be pho- 
tometered and (5) extraction of profile photometry 
from the subtracted images. 

All computations were performed with the 
frames internally subdivided into four sections 
(sub_x=sub_y=2). Differential brightness varia- 
tions of the background were fit with a second 
degree polynomial (deg_bg=2) . A convolution ker- 
nel varying quadratically with position was used 
(deg_spatial=2). The psf width (psf_width) 
was set to 33 pixels and the photometric radius 
(radphot) to 5 pixels. The reference images were 
constructed from 25 best exposures in R and 16 
in V. 

We slightly modified the reduction pipeline de- 
scribed in Paper I by introducing a procedure to 
remove photometry from epochs where a star was 
located in the proximity of bad columns. This 
task was somewhat complicated by the fact that 
the original interp program uses spline functions 
to remap each image to the template's (x, y) coor- 
dinate grid. If an image is masked before transfor- 
mation, masked regions will spread over adjacent 
columns in the remapped image. To avoid this 



problem, we performed a linear transformation of 
the bad pixel masks for each image using the coef- 
ficients output by the f itn program. The shifted 
masks were applied to subtracted images. The 
Cphot program was modified, so that it ignored 
epochs where a bad pixel was within radphot pix- 
els of a star's centroid. 

3.2. Calibration 

The transformations of instrumental magni- 
tudes to the standard system were derived from 
observations of 15, 17, 15 and 17 stars on Chips 1- 
4, respectively, in four Landolt (1992) standard 
fields, collected on 2002 May 18. Transformations 
in the following form were adopted: 

v = V + ai + a 2 (V -R) + a 3 (X - 1.25) 
v - r = h + b 2 (V -R) + b 3 (X - 1.25) 
r = R + ci + c 2 (V -R) + c 3 (X - 1.25) 

where X is the airmass. Table 1 lists the coeffi- 
cients ai, bi, Ci and the rms scatter between the 
observed and calculated standard VR magnitudes. 

These coefficients were used to calibrate the 
photometry from the images of the cluster taken 
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R [mag] R [mag] 

Fig. 3. — Left panel: The rms scatter of the i?-band light curves for stars on Chip 3 with at least 650 data 
points. The continuous curve indicates the photometric precision limit due to Poisson noise of the star and 
average sky brightness. Right panel: The detection efficiency of 0.95, 1.15, 1.30 and 1.50 Rj planets as a 
function of magnitude and rms scatter (white: 0%, black 100%), determined in §5.3. 



during the same night as the standards. The 
magnitudes from the reference images were trans- 
formed using the same color and extinction coeffi- 
cients. The offsets were determined relative to the 
calibrated photometry from cluster images taken 
on the standard night. Figure 2 shows the cal- 
ibrated V/V-R color-magnitude diagram (CMD) 
for the Chip 3 reference image 2 . 

A comparison of our T^-band magnitudes with 
the photometry of Stetson, Bruntt & Grundahl 
(2003) reveals offsets of 0.048, 0.027, 0.047 and 
0.009 mag in Chips 1-4, based on 20, 303, 3423 
and 280 stars above V = 20, respectively. We also 
find an offset of 0.022 mag in V between our Chip 
3 and Mochejska et al. (2003) data. 

The VR light curves were converted from differ- 
ential flux to instrumental magnitudes using the 
method described in Paper I. Instead of Eq. (1) 
from Paper I, we used the following relation to 
compute the total flux corresponding to the i-th 

2 The calibrated VR photometry is available from the au- 
thors via the anonymous ftp on cfa-ftp.harvard.edu, in the 
/pub/bmochejs/PISCES directory. 



image, cf. 

Ci = C re f - ACj (1) 

where Aa = c re f — a is the flux on the i-th sub- 
tracted image and c re / is the total flux on the 
reference image. This method should yield more 
accurate results because it is based on the refer- 
ence image which has a higher S /N ratio than the 
template image used previously. The instrumen- 
tal magnitudes were transformed to the standard 
system by adding offsets, computed individually 
for each star, between the instrumental and cali- 
brated reference image magnitudes. 

3.3. Astrometry 

Equatorial coordinates were determined for the 
i?-band reference image star lists. The trans- 
formation from rectangular to equatorial coordi- 
nates was derived using 964, 1012, 1476 and 951 
transformation stars from the USNO B-l catalog 
(Monet et al. 2003) in Chips 1 through 4, respec- 
tively. The mean of the absolute value of the de- 
viation between the catalog and the computed co- 
ordinates for the transformation stars was 0713 in 
right ascension and O'.'ll in declination. 
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4. Search for Transiting Planets 

4.1. Further Data Processing 

We rejected 157 i?-band epochs where less than 
5000 stars were detected on Chip 3 by DAOphot 
(Stetson 1987), and 25 other bad quality images 
from three nights. This left us with 936 highest 
quality i?-band exposures with a median seeing of 
2'.'1. We also removed 8 \^-band images, which left 
us with 318 exposures with a median seeing of 2'.'3. 

In the light curves we noticed the presence of 
offsets between different runs. This may be due 
to the periodic UV flooding of the CCD camera, 
which alters its quantum efficiency as a function of 
wavelength. To prevent the transit detection algo- 
rithm from mistaking these changes in brightness 
for transits, we added offsets between the runs, in- 
dividually for each light curve, so that the median 
magnitude was the same during each run. There 
were nine runs, each spanning from 44 to 187 data 
points. Typical sizes of the offsets were 0.008 mag 
for stars below R — 18 and 0.018 mag for stars 
between R = 18 and 19. As described in §5.6, 
this procedure greatly improves our detection effi- 
ciency. 

The left panel of Fig. 3 shows the rms scatter of 
the i?-band light curves for stars on Chip 3 with at 
least 650 data points. The continuous curve indi- 
cates the photometric precision limit due to Pois- 
son noise of the star and average sky brightness. 
The right panel shows the detection efficiency of 
0.95, 1.15, 1.30 and 1.50 Rj planets as a function 
of magnitude and rms scatter (white: 0%, black 
100%), determined in §5.3. 

4.2. Selection of Transiting Planet Candidates 

For further analysis we selected stars with at 
least 650 out of 936 good epochs and light curve 
rms below 0.05 mag. Stars above the main se- 
quence turnoff (R = 17) were rejected due to their 
large radii and, hence, very small expected transit 
depths (below 0.4%). This left us with 3074 stars 
on Chip 3, and 2975 on Chips 1, 2 and 4 (825, 
1091 and 1059 stars, respectively). 

To select transiting planet candidates we 
used the box-fitting least-squares (BLS) method 
(Kovacs, Zucker, & Mazeh 2002). Adopting a cut- 
off of 6 in Signal Detection Efficiency (SDE) and 
9 in effective signal-to- noise ratio (a), we selected 



P=2.25" R p =1.50 R, i=84.980199 - AF = 0.0139 mag 
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phase 

Fig. 4. — The original, model and combined light 
curves (upper, middle and lower panels, respec- 
tively) for a star with R = 18.01 and a planet 
with a period of 3.25 days, radius of 1.3 R.j and 
inclination of 88°. 



185 candidates on Chip 3 and 39 on Chips 1, 2 
and 4 (12, 16 and 11 candidates, respectively). 
We found three candidates on Chip 3, which were 
rejected as false detections upon closer examina- 
tion. They had similar coordinates on the image 
and their periods were all nearly exact integral 
multiples of 0.9244 days. We have found 13 other 
stars within a distance of 50 pixels whose periods 
were also such multiples. An examination of the 
period distribution of all stars revealed significant 
peaks at 4, 6 and 8 times 0.9244 days. We did not 
find any other good transiting planet candidates. 

5. Estimate of the Number of Expected 
Detections 

The number of transiting planets we should ex- 
pect to find, Np, can be derived from the following 
equation: 

N P = N.fpD (2) 

where N* is the number of stars with sufficient 
photometric precision, fp is the frequency of plan- 
ets within the investigated period range and D is 
the detection efficiency, which accounts for ran- 
dom inclinations. In §§5.1, 5.2 and 5.3 we deter- 
mine fp, A* and D. 
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Period [d] Period [d] 

Fig. 5. — Detection efficiency of transiting planets as a function of their period, relative to planets with 
all inclinations (left panel) and all transiting planets (right). Shown are the distributions for all transiting 
planets (hashed histogram), detections in the model light curves (open) and marginal (gray) and firm (solid) 
detections in the combined light curves. 



5.1. Planet Frequency 

The frequency of planets is known to increase 
with the host star's metallicity. From Figure 7 in 
Santos et al. (2004), the frequency of planets for 
stars with [Fe/H] = +0.3 - +0.4 dex is -28%, and 
it drops to ~2.5% for metallicities below [Fe/H] = 
+0.1 dex. 

The percentage of planets with periods below 
10 days is 14.6% in the Santos et al. (2004) sam- 
ple. As of 16 November 2004, the corresponding 
fractions for the planet lists on exoplanets.org 3 
and The Extrasolar Planets Encyclopaedia 4 were 
15.3% and 15.8% (excluding planets detected via 
transits). In further analysis, we adopt the value 
of 15% as the fraction of planets with periods be- 
low 10 days. 

Combining these two numbers yields /p=4.2% 
for the high metallicity cluster stars and 0.375% 
for field stars. Please note that the latter estimate 
is considerably lower than the commonly adopted 
frequency of 1%. 

3 http://exoplanets.org/planet_table.txt 
4 http://www.obspm.fr/encycl/catl.html 



5.2. The Number of Cluster and Field Stars 

Most of the cluster is contained on Chip 3 
but its main sequence (MS) is also discernible on 
Chips 1, 2 and 4. To obtain a rough estimate of 
the number of stars belonging to the cluster, we 
determined the MS fiducial line and counted as 
members all stars within 0.06 mag of it in V—R, 
on all four chips. This gives 246, 381, 2201 and 350 
"cluster" stars and 577, 710, 852 and 706 "field" 
stars on Chips 1-4, respectively. Twenty six stars 
did not have V-band data, and we assumed that 
they belong to the field. There are a total of 3178 
"cluster" and 2871 "field" stars. There are more 
"field" stars on Chip 3 than on the other chips, 
which means that some of them belong to the clus- 
ter and our color cutoff is not too liberal. On the 
other hand, a small fraction of the "cluster" stars 
are field stars, so these two biases should cancel 
out to some extent. 

5.3. Detection Efficiency 

In order to characterize our detection efficiency, 
we inserted model transits into the observed light 
curves, and tried to recover them using the BLS 
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Fig. 6. — Detection efficiency of planetary transits as a function of their inclination, relative to planets with 
all inclinations (left panel) and all transiting planets (right). 



method. 

5.3.1. Model Transit Light Curves 

The model transit light curves were defined by 
five parameters: the transit depth, AF, total tran- 
sit duration, tx, transit duration between ingress 
and egress, tp (the "flat" part of the transit), the 
period of the planet, P and the limb darkening 
coefficient, u. 

The first three parameters were computed from 
equations (1), (15) and (16) in Seager & Mallen- 
Ornelas (2003): 




Equations (4) and (5) are valid for tr^P <C 1. 
The radius of the planetary orbit, a, can be derived 
from the star's mass, M», and Kepler's third law, 
with the planet's mass M p <C M»: 



The radius and mass of the star, i?„ and M„, 
were interpolated, as a function of absolute R- 
band magnitude, Mr, from the highest metallicity 
(Z = 0.03) 7.943 Gyr isochronc of Girardi et al. 
(2000). A distance modulus (m - M) R = 13.36 
mag was used to bring the observed i?-band mag- 
nitudes to the absolute magnitude scale (Chaboyer 
ct al. 1999). 

The effects of limb darkening were simulated 
using the linear approximation first introduced by 
Milne (1921): 

I{n) = (1 - u(l - u)) (7) 

where u is the limb darkening coefficient, /i = 
cos(6), 8 is the angle between the line of sight and 
the emergent flux, and 1(1) is the intensity at the 
center of the disk. We used the grid of i?-band 
limb darkening coefficients, given as a function 
of gravity, G, and temperature, T e ff, by Claret, 
Diaz-Cordoves & Gimenez (1995). For each star, 
its G and T e // were determined from Girardi et 
al. (2000) isochrones and u was interpolated from 
four closest points in G and T e ff in the Claret et 



7 




Fig. 7. — Detection efficiency of planetary transits as a function of their radius, relative to planets with all 
inclinations (left panel) and all transiting planets (right). 



al. (1995) grid. 

In addition to P, the equations contain two 
other free parameters: the planet radius, Rp and 
the inclination of the orbit, i. A fourth parameter 
which affects the detectability of a planet is the 
epoch of the transits, To- 

5.4. Test Procedure 

We investigated the range of parameters speci- 
fied in Table 2, where P is expressed in days, Rp in 
Jupiter radii (Rj), To a s a fraction of period. We 
examined the range of periods from 1.05 to 9.85 
days and planet radii from 0.95 to 1.5 Rj, with 
a resolution of 0.2 days and 0.05 Rj, respectively. 
For To we used an increment of 5% of the period, 
and a 0.025 increment in cosi. The total number 
of combinations is 432000. 

For each combination of parameters, a random 
star was chosen without replacement from the 
sample of 3074 stars on Chip 3. When the sam- 
ple was exhausted, it was reset to the original list. 
The "observables" AF, tp and tp were computed 
and when tp > 0.5 h two light curves were gener- 
ated: the model transit light curve, and the ob- 
served light curve combined with the model (here- 
after referred to as the combined light curve). 



Figure 4 shows the original, model and combined 
light curves (upper, middle and lower panels, re- 
spectively) for a star with R = 17.36 and a planet 
with a period of 2.25 days, radius of 1.5 Rj and 
inclination of 85°. The amplitude of the transit is 
0.0139 mag, and the mass and radius of the star, 
taken from the models, are 1.03 M Q and 1.28 i? Q . 

To assess the impact of the procedure to cor- 
rect for offsets between the runs on our detection 
efficiency, we investigated three cases, where the 
correction was applied: 

A. after inserting transits, 

B. before inserting transits, 

C. was not applied at all. 

Case (B) will give us the detection efficiency if 
our data did not need to be corrected, and case 
(C) if we did not apply the corrections. Case (A) 
will give us our actual detection efficiency, and its 
comparison with cases (B) and (C) will show how 
it is affected by the applied correction procedure. 

This required us to run two sets of simulations: 
on the original (cases A and C) and corrected (case 
B) light curves. In both simulations the same list 
of parameter and star combinations was used. 
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Amplitude [mag] Amplitude [mag] 

Fig. 8. — Detection efficiency of planetary transits as a function of their amplitude, relative to planets with 
all inclinations (left panel) and all transiting planets (right). 



5.5. Detection Criteria 

An examination of the frequency of recovered 
periods, relative to the input period, revealed that 
only the peaks at 1, 2 and i P inp are distinct. 
Other aliases blend in with the background of the 
incorrectly recovered periods, so we have disre- 
garded them. 

A transit was flagged as detected if: 

1. The period recovered by BLS was within 2% 
of the input period P inp , 2 P mp or | P rnp , 

2. The BLS statistics were above the following 
thresholds: SDE > 6, a > 9. 

These detections will be referred to hereafter as 
firm. Detections where only condition (1) was ful- 
filled will be called marginal. 

5.6. Detection Efficiency 

The results of the tests are summarized in Ta- 
ble 3, which lists the test type (A-C), the number 
and percentage of transits with tx > 0.5 h (out of 
the 432000 possible parameter combinations), and 
the numbers and percentages (relative to the total 
number of transits in column 2) of transits de- 
tected in the model light curves, and of marginal 
and firm detections in the combined light curves. 



Figures 5-8 show the dependence of the detec- 
tion efficiency on period, inclination, planet radius 
and transit amplitude. The hashed, open, gray 
and solid histograms denote distributions for all 
transiting planets, planets detected in the model 
light curves, and marginal and firm detections in 
the combined light curves, respectively. Left pan- 
els show the frequency of transits and transit de- 
tections relative to planets with all inclinations. 
Right panels show the detection completeness nor- 
malized to all transiting planets (plotted as hashed 
histograms in left panels). 

The tests show that 10% of planets with peri- 
ods 1-10 days will transit their parent stars. This 
frequency drops from —24% at P = l d to —5% at 
P = 10 d . All planets with inclinations 87 - 90° 
transit their host stars, and this fraction drops to 
-80% for i = 86° and -5% for i = 78°. The 
frequency of transits increases very weakly with 
planet radius. The distribution of transit ampli- 
tudes has a wide peak stretching from 0.6% to 2%, 
centered on ~1.3%. 

The percentage of detections for the model light 
curves illustrates the limitation imposed on our de- 
tection efficiency by the temporal coverage alone. 
Due to incomplete time sampling, we are restricted 
to 75% of all planets with periods between 1 and 10 



9 



days. For periods below 4 days, our temporal cov- 
erage is sufficient to detect ~90% of all transiting 
planets, and drops to ~50% at P — 9 days. The 
detection completeness increases with decreasing 
inclination because at lower i only short period 
planets can transit their host stars. It does not 
depend on the planet radius, and it decreases with 
increasing transit amplitude. 

The source of the dependence of the detec- 
tion completeness on transit amplitude is not as 
straightforward as for the other correlations. The 
amplitude depends on the radii of the star and 
planet. Since the detection completeness was 
found to be largely independent of the planet ra- 
dius, the observed trend must stem from its depen- 
dence on the host star's radius, which is a function 
of its magnitude. Such a correlation is indeed ob- 
served, with completeness increasing for brighter 
stars (not shown here). The link between the tem- 
poral coverage and magnitude comes from the ob- 
served increase in the number of points in the light 
curve with decreasing magnitude. 

For cases A, B and C, we marginally detect 
20%, 21% and 13% of all transiting planets, and 
firmly detect 10%, 11% and 4.6%, respectively. 
Transiting planets with firm detections constitute 
83%, 84% and 64% of all stars with SDE > 6 and 
a > 9. Adding offsets between runs (case A) de- 
creases the number of firm detections by 7%, com- 
pared to the desired case, where no offsets would 
be required (case B). If the offsets were not cor- 
rected (case C), we would detect only 46% of the 
transiting planets detected in case A. 

The detection completeness for firm detections 
peaks at 20% for periods 1 — 2 d and decreases with 
period more steeply than model detections. It 
does not show a marked dependence on inclina- 
tion, and strongly increases with increasing planet 
radius, from below 2% at 1 Rj to over 20% at 
1.5 Rj. This is also apparent in the right panel 
of Fig. 3, which shows the detection efficiency of 
0.95, 1.15, 1.30 and 1.50 Rj planets as a function 
of magnitude and rms scatter (white: 0%, black 
100%). 

The detection efficiency peaks at an amplitude 
of ^1%, due to the most favorable ratio between 
the transit amplitude and photometric accuracy 
for this amplitude/magnitude range. 

The efficiency of firm transiting planet detec- 



tions, relative to planets with all orbital inclina- 
tions, L>, is 4323/432000= 1.0%. 

5.7. Number of Transiting Planets Expected 

In §§ 5.1-5.3 we determined the planet fre- 
quency fp to be 4.2% for cluster stars and 0.375% 
for field stars, the number of stars in the cluster 
and field as 3178 and 2871 and our detection effi- 
ciency D to be 1.0%. We should thus expect 1.34 
transiting planets in the cluster and an additional 
0.11 of a planet among field stars. 

5.8. Discussion 

Figure 5 demonstrates that our temporal cover- 
age is not the limiting factor. To increase the num- 
ber of expected planets it would be necessary to 
improve the photometric precision. The weather 
and seeing conditions turned out to be inferior to 
what we were expecting. A better quality CCD 
and a telescope with a larger diameter and/or bet- 
ter observing conditions would be required to im- 
prove the chances for a successful transiting planet 
search in NGC 6791. 

The estimate of 1.45 expected transiting plan- 
ets is not high enough to enable us to draw any 
conclusions from the fact that we have not de- 
tected any such events. 

The precision of this estimate is largely limited 
by the uncertainty in one of our basic assump- 
tions - the distribution of planetary radii. This 
distribution is not precisely known, and chang- 
ing it will have a marked effect on the final re- 
sult. Adopting a distribution of planetary radii 
from 1.0 to 1.35 Rj, corresponding to the radius 
range spanned by the six known transiting plan- 
ets (Konacki et al. 2004) would lower D from 1.0% 
to 0.7%. This translates to 1.08 detections, com- 
pared to 1.45 with the original radius distribution 
- a 26% decrease. 

In Paper I we made the assumption that the 
planetary radii would span the range 1-3 Rj, 
based on the radius of 1.347 Rj for the only known 
transiting planet at the time, HD 209458b (Brown 
et al. 2001). A simulation for planets in the radius 
range 1.5-3.0 Rj shows that 11% of them transit 
their parent stars, 75% are detected in the model 
light curves, and 56% and 45% are marginal and 
firm detections in the combined light curves. As- 
suming that planet radii are distributed evenly 



10 



V101 



P = 0.3348 




R V102 



P = 0.3667 



R V103 P=0.6479 



R V104 
18.5 F 1 ' 



P=1.0911 



0.5 1 

V105 P=1.2632 



20 




; 












20.2 










■■" > " 


20.4 











VI 1 1 



0.5 1 
P= 1.5487 



- 1 1 1 1 1 1 1 1 1 1 1 




















: i .... i .... i 





0.5 1 
B8 P=11.2601 



0.5 



0.5 



0.5 



0.5 1 



phase 



Fig. 9. — The i?-band light curves of nine new eclipsing binaries, four new other periodic variables, and 
cataclysmic variables B7 and B8. V107 is phased with two detected periods. 



between 1 and 3 Rj would give the percentage 
of firm detections of 37% and detection efficiency 
D = 4.0%, which translates into 5.33 expected de- 
tections in the cluster and 0.43 in the field. Our 
lack of detections does not favor such large plan- 
etary radii, in agreement with observations (Fig. 
5 in Konacki et al. 2004) and current models (Bo- 
dcnheimer, Laughlin & Lin 2003; Burrows et al. 
2004; Chabrier et al. 2004, Kornct et al. 2005). 



6. Previously Reported Candidates 

We examined the light curves of the transit- 
ing planet candidates reported by Bruntt et al. 
(2003). None of them were found to exhibit con- 
vincing periodic transits or eclipses. Even with 
our lower photometric accuracy, in most cases we 
should have detected a periodicity, if we had ob- 
served several transits. The rms of our R-band 
light curves is 0.009 for T10 and 0.003 - 0.006 for 
the remaining candidates. 
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7. Variable Stars 

We also extracted the light curves of known 
variable stars and searched for new ones by run- 
ning BLS in the period range 0.1-10 days. In Ta- 
bles 4-6 we list their revised parameters. We note, 
for the record, that variables V85, V88 and V96, 
reported as new discoveries by Bruntt et al. (2003), 
are the same as V76, V77 and V56 reported ear- 
lier by Kaluzny (2003) and in Paper I. We have 
reclassified V40 and V41 as eclipsing binaries. 

We have discovered 14 new variables: 1, 4, 7 
and 2 on Chips 1-4, respectively. Their parameters 
are listed in Tables 4 and 6 and their light curves 
are shown in Figures 9 and 10. 5 They are also plot- 
ted on the CMD in Fig. 2. Variables V101-V109 
are eclipsing binaries: V101, V102, V104-V106 are 
W UMa type contact systems, V103, V107-V109 
are detached or semi-detached binaries. In addi- 
tion to eclipses, V107 displays out of eclipse vari- 
ability with a period almost exactly twice as long 
as the orbital period. In Fig. 9 it is shown phased 
with both periods. For clarity, points in the eclipse 
have been removed from the light curve phased 
with the longer period. The other periodic vari- 
ables, V110-V113, are most likely spotted stars. 
The shape of the light curve of VI 13 varies no- 
ticeably with time. 

The non-periodic variable VI 14 is located 
slightly redward of the base of the red giant 
branch. If it belongs to the cluster, it might be a 
member of the recently proposed class of variable 
stars termed "red stragglers" (Albrow et al. 2001) 
or "sub-subgiant stars" (Mathieu et al. 2003). 
Thus far, the origin and evolutionary status of 
these stars remains unknown. 

In Fig. 10 we show the light curves of the cat- 
aclysmic variables B7 and B8. In Fig. 9 they are 
phased with cycle lengths 17.6 and 11.3 days, re- 
spectively. Equally good fits to the B7 and B8 
data are given by periods of 22.12 and 17.75 days, 
closer to the cycle lengths of 25.41 and 17.73, re- 
ported by Mochejska, Stanek & Kaluzny (2003). 
More observations are required to firmly establish 
the cycle lengths of these variables. 

In Fig. 11 we show the light curves of three 

5 The VR band photometry and finding charts for all vari- 
ables are available from the authors via the anonymous ftp 
on cfa-ftp.harvard.edu, in the /pub/bmochcjs/PISCES di- 
rectory. 



RS CVn type binaries V9, VI 2 and V33, plotted 
separately for four time intervals: 2001 July 9 - 
August 1, 2002 May 14 - June 28, 2002 Septem- 
ber 10 - November 18 and 2003 May 5 - July 10. 
The shape of their light curves varies with time, 
especially for V12 and V33. 

8. Conclusions 

In this paper we have performed an extensive 
search for transiting planets in the very old, pop- 
ulous, metal rich cluster NGC 6791. The cluster 
was monitored for over 300 hours during 84 nights. 
We have not detected any promising transiting 
planet candidates. Assuming a planet frequency 
from radial velocity surveys, we estimate that we 
should have detected 1.5 transiting planets with 
periods between 1 and 10 days, with our photo- 
metric precision and temporal coverage. The main 
limitation on our detection efficiency was imposed 
by the photometric precision. 

We have discovered 14 new variable stars in 
NGC 6791: nine eclipsing binaries, four other 
periodic variables and one non-periodic variable, 
bringing the total number of variables in this clus- 
ter to 111. We have also presented high photo- 
metric precision light curves, spanning two years, 
for all previously known variables. Many of them 
show changes in light curve shape, i.e. V9, V12 and 
V33 (Fig. 11). This phenomenon is most likely due 
to the evolution of magnetic spots on the surface 
of these stars. 

Transiting planets have proven to be more 
challenging to detect than initially expected, 
as shown by the paucity of detections from 
the many searches under way in open clusters 
(i.e. Bruntt et al. 2003; UStAPS: Street et al. 
2003; EXPLORE/OC: von Braun et al. 2004; 
STEPSS: Burke et al. 2004) and in the Galactic 
field (i.e. EXPLORE: Mallen-Ornclas et al. 2003; 
OGLE: Udalski et al. 2002a; STARE: Alonso et 
al. 2003; HAT: Bakos et al. 2004 6 ). To date, only 
six planets have been discovered independently by 
transit searches, all of them in the field, and five 
of those were initially identified by OGLE (Udal- 
ski et al. 2002a, 2002b, 2002c, 2003; Alonso et al. 



'For a more complete list of transit- 
ing planet searches, please refer to 
http: //star-www. st-and. ac .uk/~kdhl /transit s/table .html 
and http : //www . obspm . f r/encycl/searches . html 
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Table 1 
Calibration Coefficients 
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2 
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0.1664 


0.009 


0.3349 
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3 
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4 
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2.6339 
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Table 2 
Parameter Range 



Parameter 


min 


max 


step 




P (days) 


1.05 


9.85 


0.200 


45 


R P (Rj) 


0.95 


1.50 


0.050 


12 


To 
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20 


cos i 
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Table 3 

Artificial transit test statistics 



test 


all transits 


model 


marg 


inal 


firm 


type 


N 


% 


N 


% 


N 


% 


N % 


A 


43371 


10.0 


32406 


74.7 


8737 


20.1 


4323 10.0 


B 


43367 


10.0 


32380 


74.7 


9173 


21.2 


4659 10.7 


C 


43371 


10.0 


32406 


74.7 


5754 


13.3 


1998 4.6 



15 



Table 4 

Eclipsing binaries in NGC 6791 
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37 


29 


07. 


4 


19.3807 


16 


.314 


16 


.888 


0.467 





.561 
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Table 5 

Other periodic variables in NGC 6791 



ID a 200 o [h] <5 200 o [°] P [d] (R) (V) A R A v 



V42 


19 


21 


00 


2 


37 


42 


■53 


4 


o 


5064 


18 


954 


19 


597 


g 


035 


g 


.044 


v yo 


1 


21 


uu. 


2 




17 


OS 

UO 


4 


o 
u 




1 fi 

J-U 


.4/0 


I f\ 

1U 


Q9^ 
yzo 


o 

u 


00^ 
uuo 


n 
u 


004 


V110 


19 


21 


05 


g 


37 


44 


30 


4 


4 


1329 


17 


.342 


17. 


828 


o 


005 


g 


.006 


Vlll 


1 

±y 


90 
zu 


L.) 




^7 
o / 


/IS 


A ^ 
lO 


7 


4 


^AS7 


1 

±y. 


R79 
u t z 


90 
zu. 


ooo 


o 

u 


0^ 
uoo 


o 

u 


0^0 
. uou 


V84 


19 


20 


47. 


7 


37 


44 


58 


2 


4 


6258 


18 


989 


19 


836 


o 


021 


o 


027 


V44 


19 


19 


37. 


^ 


37 


41 


41 


7 


2 


2544 


17 


782 


18 


.410 


o 


013 


g 


019 


V16 


19 


21 


07. 


6 


37 


48 


09 


.6 


2 


2664 


17 


.276 


17. 


850 


o 


.019 


g 


024 


V76 


19 


20 


49. 


9 


37 


45 


50 


.9 


4 


0924 


17 


.585 


18 


270 





030 


g 


.034 


V49 


19 


20 


30 


9 


37 


36 


51 


2 


4 


9923 


15 


266 


15 


802 


g 


.007 


g 


.007 


V45 


19 


20 


16 


\ 


37 


42 


05 


9 


5 


0883 


16 


463 


17 


083 


g 


007 


g 


010 




1 

±y 


21 


1 

±y. 


O 
u 


°,7 
O 1 


1 7 


JU 


4 


r: 
o. 


1 9S7 

. ±ZO 1 


1 7 


Q^O 
. you 


1 s 


fiSS 
uoo 


o 
u 


O 4 ^ 
uoo 


o 

u 


0^7 
. UO 1 


V47 


19 


19 


39 




37 


32 


10 


g 


5 


6066 


18 


.737 


19 


978 


o 


024 


o 


.011 


V91 


19 


21 


00 


5 


37 


48 


40 


.6 


5 


6411 


17 


611 


18 


097 


g 


004 


g 


004 


VAS 
V *±o 


1 

±y 


21 


u / . 


r: 
u 


O 1 


A °, 
10 


Ofi 


« 
.u 


r: 
o. 


SOI Q 
ouiy 


1 7 


on*? 
. uuo 


1 7 


^s 


n 

u 


01 s 


o 
u 


01 ^ 

. U-LO 


V50 


19 


20 


35. 


2 


37 


31 


04 


.3 


5 


.8812 


16 


.055 


16. 


.544 





Oil 





.008 


V89 


19 


20 


56 


6 


37 


46 


36 


.2 


6 


.1577 


18 


.226 


19. 


078 





.025 





.053 


V17 


19 


20 


38. 


9 


37 


49 


04 


.5 


6 


.3656 


17 


.279 


17. 


.949 





.020 





.020 


V52 


19 


21 


20. 


9 


37 


46 


19 


.2 


6 


.9933 


17 


.016 









,005 






V51 


19 


21 


12 


2 


37 


44 


54 


.7 


7 


0315 


19 


.257 


19. 


.971 





.032 





.025 


V77 


19 


20 


52, 


9 


37 


46 


36 


.9 


7. 


1810 


16 


.206 


16. 


744 





.003 





.004 


V53 


19 


21 


00 


8 


37 


44 


35 


.4 


7. 


.1822 


18 


.257 


18. 


.803 





Oil 





.006 


V83 


19 


20 


46. 


4 


37 


44 


14 


.1 


7 


.2915 


18 


.683 


19. 


.392 





.013 





.017 


V82 


19 


20 


39. 


.7 


37 


47 


36 


.0 


7. 


.4983 


18 


.522 


19. 


.064 





.018 





.008 


V54 


19 


21 


18. 


.7 


37 


43 


36 


.4 


8. 


3141 


15 


.929 


16. 


.524 





.010 





.010 


V98 


19 


20 


56. 


5 


37 


45 


38 


.7 


8. 


.3405 


16 


.418 









.003 






V112 


19 


20 


04. 


2 


37 


48 


33 


.4 


8. 


.7902 


14 


.947 


15, 


.467 





.003 





.002 


V95 


19 


20 


43. 


1 


37 


47 


32 


.5 


9 


.6832 


18 


.509 


19. 


147 





.014 





.010 


V65 


19 


20 


52. 


5 


37 


47 


30 


.5 


11 


1091 


15 


.645 


16. 


.272 





.004 





.006 


V56 


19 


20 


45. 


3 


37 


45 


48 


.8 


12 


.3832 


16 


.518 


17 


.081 





.006 





.003 


V57 


19 


20 


57. 


9 


37 


31 


07 


.0 


13 


1536 


17 


.460 


18. 


.359 





.010 





.007 


V58 


19 


21 


14. 


.5 


37 


48 


04 


.4 


13 


.2597 


17 


.028 


17. 


.544 





.015 





.016 


V113 


19 


20 


34. 


9 


37 


48 


14 


.8 


13 


.5908 


16 


.995 


17. 


.568 





.009 





.013 


V97 


19 


20 


49. 


2 


37 


49 


14 


.8 


13 


6206 


15 


.813 


16 


.517 





.003 





.004 


V59 


19 


20 


21 


5 


37 


48 


21 


.9 


13 


8331 


17 


.221 


17 


.781 





.056 





.049 


V64 


19 


21 


11. 


4 


37 


29 


55 


.4 


14 


.2427 


15 


.882 


16. 


.473 





.005 





.012 


V81 


19 


20 


49. 


.7 


37 


48 


08 


.7 


16 


.6182 


16 


.373 


16 


.905 





.003 





.005 


V100 


19 


21 


01. 


.8 


37 


45 


41 


.9 


23 


.9468 


16 


.544 


17, 


.151 





.024 





.021 


V66 


19 


21 


08. 


4 


37 


44 


55 


.2 


50. 


.4976 


15 


.372 


16. 


.119 





.068 





.074 


V71 


19 


21 


10. 


5 


37 


43 


24 


.8 


51 


.9430 


16 


.563 


17 


.291 





.057 





.067 


V67 


19 


21 


03. 


.7 


37 


48 


03 


.7 


66. 


.7944 


16 


.180 


16 


.959 





.056 





.071 
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Table 6 

Miscellaneous variables in NGC 6791 



ID 


02000 [h] 


feooo [°] 


Rmax 


Vmai 


A R 


A v 


B07 


19 21 07.4 


37 47 56.5 


17.448 


17.581 


0.819 


0.919 


B08 


m on qc; *7 

19 2U oo.Y 


37 44 52.3 


18.820 


18.716 


1.896 


3.403 


V10 


19 21 11.8 


37 47 58.1 


18.922 


19.642 


0.103 


0.435 


V21 


19 20 57.3 


37 45 36.9 


16.997 


17.547 


0.016 


0.041 


V62 


19 21 03.0 


37 43 51.8 


18.618 


19.189 


0.105 


0.214 


V63 


19 19 40.0 


37 29 45.1 


16.317 


17.058 


0.024 


0.066 


V70 


19 20 32.2 


37 44 21.0 


99.999 


14.722 


0.000 


0.427 


V74 


19 21 07.2 


37 44 34.9 


99.999 


14.660 


0.000 


0.025 


V75 


19 20 47.9 


37 45 58.8 


16.829 


17.374 


0.021 


0.053 


V79 


19 20 55.2 


37 46 39.7 


18.034 


18.631 


0.083 


0.261 


V86 


19 20 50.1 


37 48 31.7 


18.846 


19.478 


0.134 


0.550 


V87 


19 20 52.8 


37 44 58.8 


17.637 


18.189 


0.038 


0.094 


V90 


19 20 58.9 


37 44 47.1 


17.608 


18.159 


0.039 


0.148 


V94 


19 20 42.5 


37 44 36.9 


17.023 


17.563 


0.034 


0.081 


V99 


19 20 57.1 


37 48 12.1 


16.866 


17.527 


0.025 


0.045 


V114 


19 20 00.0 


37 48 44.7 


16.882 


17.601 


0.034 


0.057 
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